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Abstract. The Hamiltonian approach to the theory of dual isomonodromic defor- 
mations is developed within the framework of rational classical A-matrix structures 
on loop algebras. Particular solutions to the isomonodromic deformation equations 
appearing in the computation of correlation functions in integrable quantum field the- 
ory models are constructed through the Riemann-Hilbert problem method. The corre- 
sponding r-functions are shown to be given by the Fredholm determinant of a special 
class of integral operators. 
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1 Monodromy Preserving Hamiltonian Systems 
1.1 Isomonodromic Deformation Equations 

Monodromy preserving deformations of rational covariant derivative operators 
of the form: 

^a:=^t-AA(A), (1.1) 



where 



<9A 

" N- 

m):=B + J2^ (1-2) 

i=l 

B is the diagonal r x r matrix 

B = diag(/3 1 ,...,/3 r ), (1.3) 

and the matrices {iVj}j = i... n are r x r matrix functions of the n + r deforma- 
tion parameters {ai,/3 a }, were studied by Jimbo et. al. in (Jimbo et al (1980)), 
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(Jimbo et al (1981)), (Jimbo and Miwa (1981)). It was shown there that the most 
general differentable monodromy preserving deformations of such operators are 
determined by the intcgrablc Pfaffian system: 

n 

dNi = - Y^Wi, N^dlogiai - a,) - [N u d( ai B) + O], (1.4) 

3 = 1 

where is the r x r matrix with elements 

n 

O ab - (1 - <U)(E Ni) ab dlog(f3 a - f3 b ). (1.5) 

Such operators and their monodromy are of great importance in the theory 
of quantum integrable systems, since the computation of correlation functions 
in such systems very often leads to particular solutions to such systems (Its et 
al. (1990)), (Korepin et al. (1993)). In the following subsection, it will be shown 
how these equations may be understood as a compatible set of nonautonomous 
Hamiltonian systems generated by commuting Hamiltonians that are spectral 
invariants of the matrix W(A). In subsequent sections, the classical i?-matrix 
approach to such systems will be explained and the computation of certain so- 
lutions related to Fredholm determinant calculations via the matrix Riemann- 
Hilbcrt problem will be decribed. 

1.2 Nonautonomous Hamiltonian Structure 

We begin with the Lie Poisson structure on (©™ = i0l(r))*, defined by the following 
Poisson brackets between the various matrix elements 

{(A^U, (NjU} = S t3 (5 bc (N t ) ad - S ad (N t ) cb ) . (1.6) 

The system (1.4) can be expressed in multi-Hamiltonian form by introducing 
the following Hamiltonian 1-form on the parameter space 

n r 

i—l a—1 

where 

gi : = tr(BjV,)+V ^^) ; » = !,...,„ (1.8) 

K a := 2^ a t (Ni) aa + ^ _ J , a=l,...,r. (1.9) 

i=i t=i * 3a P b 

Equations (1.4) may then equivalently be written in multi-Hamiltonian form as 

dNi = {Ni,6}. (1.10) 
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Involutivcness of the H^s and K^s, which will be explained in the following 
section, then implies that the differential form 8 on the parameter space is in 
fact closed 

d0 = 0, (1.11) 

which allows one to introduce the r-function (cf. Jimbo et al (1980), Jimbo et 
al (1981)) by the formula 

= eflogr. (1.12) 

In the next section, we show how the above Hamiltonian structure may very 
naturally be viewed as the restriction of the rational i?-matrix structure on a 
loop algebra to a finite dimensional Poisson submanifold. From this the commu- 
tativity of the Hamiltonians (Hi, K a ) defined above follows. 

2 Loop Algebra Moment Maps, Spectral Invariants and 
Isomonodromic Deformation Equations 

The following discussion is based on the approach developed in Harnad (1994). 

2.1 Dual Moment Maps and Split ii-matrix Structure 

We introduce an auxiliary symplectic vector space (M, lu), which will be referred 
to as the generalized Moser space, consisting of pairs (F, G) of complex N x r 
matrices whose elements are viewed as canonically conjugate variables. Thus, 
the symplectic form lo is just 

oj :=tr(dF AdG T ). (2.1) 

We use the following notation to denote the loop algebra of r x r matrices 
depending on a loop parameter A, viewed as a point on a circle S 1 in the complex 
A-plane, and its splitting into negative a positive Fourier components 

fll(r)=fl[(r)++fl[(r)_~fl[(r)* (2.2) 
Ql(r)* ± ~ fl I(r) T (2.3) 

The identification with the dual space indicated in eq. (2.2) is defined through 
the Ad-invariant scalar product 

< X X ,X 2 >= <f tr (X 1 (\)X 2 (X)) dX. (2.4) 
Js 1 

The rational i?-matrix structure is obtained by just redefining the Lie algebra 
structure in such a way that the new algebra splits into a Lie algebraic direct 
sum of the positive and negative frequency parts, with a change of sign in the 
Lie product for the second summand 



5 l R (r) := fl [(r)+ fll(r)_ - fll fl (r)*. 



(2.5) 
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The rational i?-matrix structure is then just the corresponding Lie Poisson struc- 
ture on gl R (r)* . Expressed in terms of the individual matrix elements this gives 

EAUA), "M] = (2.6) 

A — fl 

This can be expressed more succintly in the tensorial (St. Petersburg) notation 
as follows 

{JV(A)®Jvv»} = [r(A - fx), W(A) ® I + I ® Af{n)] , (2.7) 

where Af(A) and A/"(/i) are viewed as endomorphisms of C" and the rational 
i?-matrix r(A — /x) is the endomorphism of C™ <g> C™ defined by 

r(A-/i) := G End(C"®C"), (2.8) 
A — /i 

where P12 denotes the endomorphism that interchanges the first and second 
factors in C" <8> C™. 

Now let A and B be the diagonal N x N and r x r matrices, respectively, 
defined by 

A := diag(aj) e gl(N), B := diag(/3 a ) £ g[(r), (2.9) 

where the eigenvalues {ctij}j = i...„ have multiplicities {fcj}j = i ; ...„, and the eigen- 
values {/3 }a=i— r are multiplicity free. Define the Poisson subspace g^ C flt fl (r)* 
by 

n N- " 
g^:={^(A)=S + ^--^-, ^e [(r)}~]T f(r). (2.10) 

Then the following defines a Poisson quotient map of the symplcctic space M, 
such that the image is identified with a Poisson submanifold of gg. 

j£ : M — > g^ 

: (F, G) 1 — ► B + G T (A - XI^F =: 7V(A) (2.11) 
" AT 

JV(A) = B + ^-^egl (r) fl (2.12) 



A — cti 



Ni := —GfFi, F l ,G l eM^ xr . (2.13) 

We may now apply the standard classical i?-matrix theory to deduce a set 
of commuting Hamiltonian^ flows on the space M, generated by the spectral 
invariant Hamiltonians on gl(r)* , pulled back to M via the above Poisson map, 
for which the Hamiltonian flows are represented by Lax equations. The resulting 
Hamiltonian flows are therefore isospectral for the matrix JV(X). We denote by 

T£:=T(0[(r)*)| g . (2.14) 

the ring of spectral invariants restricted to g^. Then the classical i?-matrix 
theory in this case tells us that: 
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(i) lg is Poisson commutative. 

(ii) For H gig, Hamilton's equations have the Lax form: 

^ = H?, A/], (2.15) 

A* := <TdH+{N) + (<j- l)dH_{M) =: V a (dH(JV)), (2.16) 

where the subscripts ± denote projections to the positive and negative Fourier 
components and a £ R is arbitrary. 

It follows that the spectral curve defined by the characteristic equation 

det(7V(A) - zl r ) = 0, (2.17) 

where 

JV(X)B + G T (A- XI N )-\ (2.18) 

is invariant under the resulting Hamiltonian flows. 

In order to apply this to the isomonodromic deformation equations considered 
above, we must adapt these results to the case of nonautonomous Hamiltonian 
systems, in which the flow parameters are reinterpreted as deformation parame- 
ters upon which the spectral invariant Hamiltonians may explicitly depend. This 
will be done in the next subsection. 



2.2 Nonautonomous systems: isomonodromic deformations 

Letting the matrices A and B depend explicitly on some deformation parameter 
t 

A = A(t), B = B(t), (2.19) 

the above Lax equations must be modified to take the resulting explicit t 
dependence of the matrix 7V(A) into account. This gives the nonautonomous 
system 

+ <»•»> 

Suppose now that, for some a and H, the following special condition holds: 

^L = ^. (2.21) 

It follows that Hamilton's equations become isomonodromic deformation equa- 
tions, since eq. (2.20) then takes the form of commutativity conditions 

[D x , V t \ = 0, (2.22) 

where 

V x :=-^-M(X) (2.23) 

V t :=^ t -A». (2.24) 
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These are precisely the necessary conditions for the invariance of the monodromy 
of the operator V x under deformations in the parameter t. 

To apply this to the system (1.4), we choose the following set of spectral 
invariant Hamiltonians {Hi £ 1g}i=\,... n 

Hi(Af) :=-?-(( tr((AA(A)) 2 dA = tr(BNi) + V t l!3M . ( 2 .25) 

jlti 

The autonomous form of Hamilton's equations that result are then 

^ = -[(dHiUAf], (2.26) 

where 

AT- ~ 

{dHi)_ = - — l — G fl[(r)_. (2.27) 

A - a,i 

Identifying the various deformation parameters now with the eigenvalues of the 
matrix A, {ti — ai}i=i t ..., n , the Lax equations are modified to the following 
form: 

These are just the commutativity conditions 

[V x ,Vi] = 0, i = l,...,n, (2.29) 
V X :=^L-M(\) (2.30) 

Vi ■■=■£- + idHt). = ■£- + t— i-. (2.31) 

don octi A — cti 

guaranteeing the preservation of the monodromy of the operator V x under the 
deformations generated by varying the a^s. Evaluating residues at {ai}i=i,...n 
gives 

^ = JM, 

^-P + t^."d- ("») 

which are just the on components of the differential system (1.4). 
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2.3 Dual Isomonodromic System 

To obtain the (3 a components of this system, it is convenient to introduce an- 
other representation, in terms of a second system of rational covariant derivative 
operators whose monodromy will also be preserved: the dual isomonodromic sys- 
tem. Define another loop algebra gl(N), consisting of N x TV matrices depending 
similarly on a loop parameter z that lies on a circle S 1 in the complex z-plane, 
with corresponding splitting into positive and negative Fourier components: 

Ql(N) = sl(N)+ + fl[(A0- ~ &KN)*, (2.34) 

and corresponding rational i?-matrix structure gl(N)* R on the dual space. We 
also define a corresponding Poisson subspace g^ C Q[{N)* R consisting of rational 
elements M(z) of the form 

gf := {M(z) = -A + J2 ^7"} ~ E A r W- ( 2 ' 35 ) 
Introduce the "dual" Poisson map from M to g^ as: 

jf : (F, G) i — ► -A - F(5 - zI N y x (F := X(z), (2.36) 

and denote the dual ring of spectral invariants restricted to g^ as 

Zf :=T(0[(AO*)| g B. (2.37) 

We then have the remarkable fact that the spectral rings 2/f and T R coincide 
when pulled back under the respective Poisson maps Jjj* and Jg to M. 

Theorem 2.1 (duality theorem) : The two spectral invariant rings J%* (2%) 
and Jg*{Ig) coincide. 

Proof: This essentially follows from the simple linear algebra identity 

det(A - XI N )dct(B + G T (A - Mn^F - zl r ) 
= det(B - zI r )det(A + F(B - zl r y x G T - XI N ), (2.38) 

which implies that the spectral curves of ftf(X) and M{z) are identical. 

Now define the set of spectral invariant Hamiltonians {K a G Xf } =i,...n 
similarly to the H^s, as: 

K a ■= i tv(M(z)) 2 dz. (2.39) 
4tti J z=0a 

On g^, these similarly generate the equations: 
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where 



(dK a )_(z) = -^-£ 3 l(N)_, (2.41) 

Z Pa 



which imply the invariance of the monodromy of the rational covariant derivative 
operator 

V z := — - M(z) (2.42) 

under the deformations generated by changes in the parameters {[3 a }a=i, - r- 

Using Theorem 2.1, we may also pull back the K a } s under the map Jg, 
to determine corresponding spectral invariant Hamiltonian functions of W(A). 
These again generate isomonodromic deformation equations for the operator V\, 
which may be expressed 

g-^M + Z&k, (2.43, 



where 

,, r/ > / \ \ > j-, . •^•^E a N i E b + E h NiE a 
(dK a )+(X) = XE a + 2^ 2^ , (2.44) 

6=1 j=l Pa > Db 

and E a denotes the elementary r x r matrix with diagonal entry 1 in the aa 
position and zero elsewhere. Evaluating residues at A = ai gives 

^ = [(dK a ) + (a i ),N i ], (2.45) 

which are precisely the /3 a components of the isomonodromic system (1.4). Sim- 
ilarly, viewing {i?i}i=i,...„ as Hamiltonians defined on g^, these generate the 
dual equations, which imply the invariance of the monodromy of the operator 

where 

«IH M *) = -*E I + ±± E < M - E < +E ' M ° E > , ,2.47) 

J = i a=l ■> 

and Ei denotes the elementary n x n matrix with entry 1 in the ii position. 
Evaluating residues at z — (i a gives these in the form 

^ = [(dH i ) + ((3 a ),M a }. (2.48) 
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3 Isomonodromic Deformations and the 
Riemann — Hilbert Problem 

In this last section, we discuss certain specific solutions of the above isomon- 
odromic deformation equations which can be constructed through application of 
the Zakharov-Shabat "dressing" method (Novikov et al. (1984)). This class of 
solutions is of particular interest from the viewpoint of applications, since they 
arise in the calculation of correlation functions for quantum integrablc systems 
(Korepin et al. (1993)), (Harnad and Its (1997)) and spectral distributions in the 
theory of random matrices (Tracy and Widom (1994)), (Harnad et al. (1994)). 
The results quoted in this section are based on joint work with A. Its; the full 
details may be found in the joint paper (Harnad and Its (1997)). 

The particular class of solutions to the isomonodromic deformation equations 
in question may be constructed by applying the dressing method, based on the 
matrix Riemann-Hilbcrt problem, suitably adapted to this case. To do this, we 
first introduce a vacuum solution <Jo, which is chosen as the invertible r x r 
matrix function obtained by exponentiating XB 

#o(A) := e XB . (3.1) 

We then introduce a family of loop group elements Hq{\), viewed as GZ (re- 
valued functions defined along some oriented, closed curve r chosen, in this 
case, to pass through the points consecutively, with the latter ordered by 

their subscript labels. We also assume in the following that the number n of 
such points is even, and write n — 2m. (If the number happens to be odd, we 
just increase it by adding A = oo as the last point.) Let {Fj}j =1 ... m denote the 
segment of r between «2j-i and ot^j and let 6j(X) denote the characteristic 
function, along r, of the interval Jj. We define H (\) as the piecewise constant 
element of the form 

n 

H (X) := I r + 2™^f,gJ^(A), (3.2) 

where {fj, Ej}j=i---m i s anv fixed set of r xp complex, rectangular matrices, with 
p < r, satisfying the null conditions 

gjf fe =0, Vj,k. (3.3) 

The relevant matrix Riemann-Hilbert problem consists of finding a nonsingular 
r x r matrix valued function x{X) that is analytic on the complement of r, 
extending to A = oo off r, with asymptotic form 

X (A) ~ I r + (A- 1 ) (3.4) 

for A — > oo, and has cut discontinuities across r given by 



X_(A) = x+ (A)ff(A), Aef, 



(3.5) 
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where x+(A) and X-(A) are the limiting values of x(A) as T is approached from 
the left and the right, respectively, and if (A) is the rxr invertible matrix valued 
function along T defined by 

H(X)=%(X)H (X)^ 1 (X). (3.6) 

Following the Zakharov-Shabat dressing method, we define the dressed wave 
function as 

*±(A) := x±(A)*b(A), (3.7) 
with limiting values <?± on either side of the segments of T given by 

*±(A) := X±(A)#b(A). (3.8) 

We then have the following result, which is quoted here from Harnad and Its 
(1997), 

Theorem 3.1: The wave function defined by (3.8) satisfies the equations 



SA 




with the Nj 's given by 
where 



& = 0, (3.9) 

, , # = 0, (3.10) 

OCtj A — Qj 

Nj := -F.-Gj, (3.11) 



Fj := lim F(A) Gj := {-l) j Urn G(A). (3.12) 

A — >ctj A — >aj 

This implies the commutativity 

[V x ,V ai ]=0, [D ai ,V aj }=0, i,j = l,...,n (3.13) 
of the operators 

P Q , :=^- + T ^, (3.15) 

and hence the invariance of the monodromy data of the operator V\ under 
changes in the parameters {ctj} . 

Thus, the operators defined in eqs. (3.14)-(3.15) represent a solution to the otj 
components of the isomonodromic system (1.4). The following result, also quoted 
from Harnad and Its (1997), shows that the same construction also provides a 
solution to the (3 a components. 
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Theorem 3.2: The wave function tf'(A) also satisfies the equations 

Vp a $ = Q, a = l,...,r, (3.16) 
where the operators {£>g }6=i,...,r are defined by 

^-ss-^-s:-* h-ir — • (3 - 17) 

Ma 

with the Ni's given by eqs. (3.11), (3.12). This implies the commutativity con- 
ditions 

[V x , Z>/j a ]=0, [D^, Vp b }=0, o,6=l,...,r, (3.18) 

and hence the invariance of the monodromy data of V\ under the deformations 
parameterized by {/3 }o=i,...,r- 

From the viewpoint of applications to quantum integrable systems (Its et 
al. (1990), Korepin et al. (1993)) and the spectral theory of random matrices 
(Tracy and Widom (1994), Harnad et al. (1994)) this construction has particu- 
lar importance, since the underlying T-function, as defined in eq. (1.12), is just 
the Fredholm determinant of an integral operator that may be constructed from 
the same data, and which gives the correlation functions and spectral distribu- 
tion generating functions in question. This result is contained in the following 
theorem, also quoted from Harnad and Its (1997). 

Theorem 3.3: Let K be the p x p matrix Fredholm integral operator acting 
on C p -valued functions v(A), 

K(v)(A) = J K(X, /i)v(/x)d/i, (3.19) 

defined along the curve T, with integral kernel given by 

K(X,») = qMEl, ( 3. 2 0) 

A — fl 

where f , g are the rectangular r x p matrix valued functions 

m 

f(A) :=<MA)]Tf^-(A) (3.21) 

3=1 

m 

g(A):=«(A))- 1 ^g/ J (A). (3.22) 

3=1 

Then the logarithmic derivative of the Fredholm determinant is given by 

n r 

dlndet(I-K) = w = ^H k da k + ^K a dp a , (3.23) 

fe=l a=l 
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where the individual factors may be expressed 



mndetq-K) £ ^ 

dlndet(I-K) 



K n = 



= E «i M)« + E — ^ • ( 3 - 25 ) 

, = 1 6=1 Pa Ph 

Hence, det(I — K) may be idcntihcd as the r-function defined in eq. (1-12). 

Finally, it should be mentioned that the dual isomonodromic systems defined 
in eqs. (2.42)-(2.48) may be derived in exactly the same way, by interchanging 
the roles of the matrices A and^-B when defining the vacuum wave function 
(3.1). The corresponding curve r must be chosen, in the complex z-plane, so 
as to pass through the parameters {/3 a } giving the diagonal elements of the 
matrix B. The resulting r-function turns out to just be given by the Frcdholm 
determinant of the dual Frcdholm integral operator K, which is related to the 
operator K appearing in Theorem 3.2 by taking a Fourier-Laplace transform 
along the curves r and f. Full details regarding this result, as well as a number 
of related results, including generalizations to isomonodromic deformations of 
operators having higher order pole singularities, may be found in Harnad and 
Its (1997). 
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